Self-consistent calculation of the electron distribution near a Quantum-Point Contact 

in the integer Quantum Hall Effect 
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In this work we implement the self-consistent Thomas-Fermi-Poisson approach to a homogeneous 
two dimensional electron system (2DES). We compute the electrostatic potential produced inside a 
semiconductor structure by a quantum-point-contact (QPC) placed at the surface of the semicon- 
ductor and biased with appropriate voltages. The model is based on a semi-analytical solution of 
the Laplace equation. Starting from the calculated confining potential, the self-consistent (screened) 
potential and the electron densities are calculated for finite temperature and magnetic field. We ob- 
serve that there are mainly three characteristic rearrangements of the incompressible "edge" states, 
which will determine the current distribution near a QPC. 

PACS numbers: 73.20.-r, 73.50.Jt, 71.70.Di 



I. 



INTRODUCTION 



A quantum point contact (QPC) is constructed by 
geometric or electrostatic confinement of a two dimen- 
sional electron system (2DES). The conductance through 
them is quantized^^ and they play a crucial role in the 
field of mesoscopic quantum transport. Their properties 
have been investigated in a wide variety of experiments, 
which include the observation of the 0.7 anomalji^, 
quantum dots coupled to QPCs^, Quantum-Hall effect 
(QHE) based Mach-ZenderiSii and Aharonov-Bohm in- 
terferometers. This has lead to extensive investigations 
of the electrostatic and transport properties of QPCs, 
both with and without a quantizing magnetic field. Many 
different techniques have been used to find the elec- 
tronic density distribution near a QPC, ranging from nu- 
merical Poisson-Schrodinger solutions^ to spin-density- 
functional theorji^ and phenomenological approachesiS. 
It has been possible to treat realistic samples mostly 
only within simplified electrostatic calculations, neglect- 
ing screening effects. On the other hand, when including 
interactions the calculations become more complicated, 
thus one usually sacrifices handling realistic geometries. 

Recent experiments have succeeded in developing and 
analyzing a QHE based electronic Mach-Zcnder inter- 
ferometer (MZI)^, making use of the integer QHE edge 
stateai as single-channel chiral quantum wires. A key 
element of these experiments are the QPCs, which play 
the role of the beam splitters of the optical setup. The 
electrostatic potential and electronic density distribu- 
tions in and near the QPCs play an important role in 
understanding the rearrangement of the edge states in- 
volved. Moreover, the electron-electron interaction has 
been proposed^ as one of the origins of dephasing in such 
an electronic MZI, such that a self-consistent calculation 
of the electrostatic potential may also be viewed as a first 
step towards a quantitative understanding of this issue. 
So far, the theoretical description of dephasing in the 
electronic MZI via classicaliiiiSiiii^ or quantum noise 
fieldsi^ii^ and other approachesii has focused on features 



supposed to be independent of its specific realization (see 
Ref. ^3 foi' * recent review). However, a more detailed 
analysis of the QHE related physics, taking account of 
interaction effects, will certainly be needed for a direct 
comparison with experimental data. In this paper, we 
will provide a detailed numerical analysis of the electro- 
statics of QPCs in the integer QHE, assuming geometries 
adapted to those used in the MZI experiment. Our work 
will produce the electron density and electrostatic poten- 
tial, based on the self-consistent Thomas-Fermi-Poisson 
approximation, to which we refer as TEA in the following. 

We would like to point out the following observation 
regarding the Mach-Zchnder experiment, where a yet- 
unexplained beating pattern observed in the visibility 
(interference contrast) as a function of bias voltage was 
surprisingly found to have a period independent of the 
length of the interferometer arms. Such a result would 
seem less surprising if all the relevant interaction physics 
leading to the beating pattern were actually taking place 
in the vicinity of the QPC. This provides strong encour- 
agement for future more detailed work on the coherent 
transport properties of these QPCs. 

Although it has been more than two decades since the 
discovery of the quantized Hall effectiS, the microscopic 
picture of current distribution in the sample and the in- 
terplay of the current distribution with the Hall plateaus 
is still under debate. In recent experiments, the Hall 
potential distribution and the local electronic compress- 
ibility have been investigated in a Hall bar geometry by 
a low-temperature scanning force microscopeSS and by 
a single-electron-transistoriSi, respectively. This has mo- 
tivated theoretical^^ work, where a self-consistent TEA 
calculation has been used to obtain electrostatic quanti- 
ties. 

Self-consistent screening calculations show that the 
2DES contains two different kinds of regions, namely the 
quasi-metallic compressible and quasi-insulating incom- 
pressible regions22i24. The electron distribution within 
the Hall bar depends on the "pinning" of the Fermi level 
to highly degenerate Landau levels. Wherever the Fermi 
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level lies within a Landau level with its high density of 
states (DOS), the system is known to be compressible 
(leading to screening and correspondingly to a flat po- 
tential profile), otherwise it is incompressible, with a con- 
stant electron density and, in general, a spatially vary- 
ing potential due to the absence of screening. Moreover, 
based on these results for the potential and density dis- 
tributions, one may employ a local version of Ohm's law 
(together with Maxwell's equations and an appropriate 
model for the conductivity tensor) to calculate the cur- 
rent distribution, imposing a given overall external cur- 
rent for the in-plane geometry . These results are mostly 
consistent with experiments except that within the self- 
consistent TFA one obtains an incompressible strip (IS) 
for a large interval of magnetic field values which leads 
to coexistence of several IS's with different local filling 
factors. Recently, this theory has been improved in two 
aspects2S*2&: (i) the finite extent of the wave functions 
was taken into account in obtaining electrostatic quan- 
tities (rather than using delta functions), (ii) the find- 
ings of the full Hartree calculations were simulated by 
a simple averaging of the local conductivities over the 
Fermi wave length , thereby relaxing the strict locality of 
Ohm's law for realistic sample sizes. A very important 
outcome of this model is that there can exist only one 
incompressible edge state at one side of the sample for a 
given magnetic field value. Indeed this is differing dras- 
tically from the Chklovskii-Shklovskii-Glazman (CSG)^ 
and the Landauer-Biittikei— picture, where more than 
one edge state can exist and is necessary to "explain" 
the QHE. In the CSG scheme a non-self-consistent TFA 
(which is called the "electrostatic approximation") was 
used. However, it is clear that if the widths of the IS's 
(where the potential variation is observed) become com- 
parable with the magnetic length, the TFA is not valid, 
thus the results obtained within this model are not reli- 
able any more. In principle, similar results to Ref. p5l | 
were reported by T. Suzuki and T. Ando2^, quite some 
time ago and recently by S. Ihnatsenka and I. V. Zo- 
zoulenko^ in the context of spin-density-functional the- 
ory. With the improvements on the self-consistent TFA 
mentioned above, together with taking into account the 
disorder potential^ and using the self-consistent Born 
approximation^ to calculate the local conductivity ten- 
sor, one obtains well developed Hall plateaus, with the 
longitudinal resistivity vanishing to a very high accu- 
racy, and one is also able to represent correctly the inter- 
plateau transition regions. Wherever one observes an IS, 
the longitudinal conductivity becomes zero, and as a con- 
sequence also the corresponding local resistance (and the 
total resistance) vanishes. Thus, according to Ohm's law, 
the current flows through the incompressible region. In 
addition, the Hall conductance becomes equal to the local 
value of the quantized conductance. Finally all the three 
experimentally observed^ qualitatively different regimes 
of how the Hall potential drops across the sample have 
been reproduced theoretically without artifacts of the 
TFA2S. The model described above has also been suc- 



cessfully applied to an electron-electron bilayer system^^ 
and provided a qualitative explanation^i of the magneto- 
resistance hysteresis that has been reported recentlji^Si^. 
For all of these reasons, we feel confident in applying this 
theory to our analysis of the MZI setup. 

Motivated by the experimental and theoretical findings 
ascertaining the importance of the interaction effects in 
the integer Quantum Hall regime, in this work we will 
show that the mutual Coulomb interaction between the 
electrons leads to interesting non-linear phenomena in 
the potential and electron distribution in close proximity 
of the QPCs. Based on the self-consistent Thomas-Fermi- 
Poisson approximation (TFA), we will consider realistic 
QPC geometries and examine the distribution of the in- 
compressible regions depending on the field strength and 
sample parameters. 

The rest of this paper is organized as follows: In SecHII 
the electrostatic potential produced by an arbitrary sur- 
face gate will be discussed, by solving the Laplace equa- 
tion without screening effects. In Sec. IIIII we review the 
TFA in a 2DES. In Sec.|lVlwe will first present the well 
known general results of the TFA for a homogeneous 
2DES at zero magnetic field B and zero temperature, 
and we will investigate the electron density and electro- 
static potential profiles of a (i) simple square gate geom- 
etry and (ii) a generic QPC, before (iii) systematically 
investigating the positions of the incompressible strips 
depending on magnetic field and geometric parameters. 
We conclude with a discussion in Sec. 



II. ELECTROSTATICS OF THE GATES 

As mentioned in the introduction, there is a tradeoff 
between simulating realistic QPC geometries and includ- 
ing the interaction effects within a reasonable approxima- 
tion. In this paper, we present an intermediate approach, 
which considers realistic QPC structures but interactions 
of the electrons are handled within a Thomas-Fermi ap- 
proximation (TFA), which is valid for relatively "large" 
QPCs (> 100 nm). One can obtain, in a semi-analytical 
fashion, the electrostatic potential generated by an arbi- 
trary metallic gate at the surface by solving the Laplace 
equation for the given boundary conditions. Afterwards, 
it is possible to obtain the electron and potential distri- 
bution in the 2DES, within the TFA, both for vanishing 
and finite magnetic fields (B), and at low temperatures 
at B > 0. 

Here we briefly summarize the semi-analytical model 
developed by J.Davies and co-workers^. The aim of this 
section is to calculate the electrostatic potential on a 
plane at some position z below the surface of the semi- 
conductor, which is partially covered by a patterned gate. 
The surface occupies the 2 = plane and z is measured 
into the material. The un-patterned surface is taken to 
be pinned so we can set the potential Vup{r, 0) = there, 
with Vgate (r,0) = Vg on the gate. We use lower-case 
letters like r = (x, y) to denote two-dimensional vectors 
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FIG. 1: The image of the QPC (gray scale). The polygons 
are used to define the gates on the two-dimensional mesh. 
The 2DES resides under the dark (grey) regions, with a bulk 
electron density of 1.7 x 10~^^cm~^ (see Ref. due to mem- 
ory limitations the quality of the figure is reduced. The white 
line scales 200nm). 



with the corresponding upper-case letters for three di- 
mensional vectors like R = {x,y,z) = (r, z). Thus the 
problem is to find a solution, VextC^), to the Laplace 
equation V^Veit = 0, given the value on the plane 
z ^ 0, and subject to the further boundary condition 
dVext/dz — > as z — > oo. One route is to start by mak- 
ing a two-dimensional Fourier transform from Vexti^i^) 
to Vext (q, 0) . The dependence on z is a decaying exponen- 
tial to satisfy Laplaces equation and the boundary con- 
dition at z = oo: Vext{(\, z) = Vext{(\, 0) exp(-|qz|). This 
multiplication of the Fourier transform is equivalent to a 
convolution in real space. Taking the two-dimensional in- 
verse Fourier transform of cxp(— |gz|) leads to the general 
result. 

14xt(r,z) = i / l"^' V.ateiy' ,Q)dv' , 

Kj 27r(z^ + |r — r'l"')-^/"' 

(1) 

where k is the dielectric constant of the considered 
hctcro-structure. Now one can evaluate the potential 
in the plane of the 2DES, z = d, for a given gate and 
potential distribution on the surface. The derivation of 
some important shapes like triangle, rectangle and poly- 
gons is provided in the work cited above, which has been 
successfully applied to quantum dot systems^. For our 
geometry, we will use the result for the polygons. 



III. ELECTRON-ELECTRON INTERACTION: 
THOMAS-FERMI-POISSON APPROXIMATION 

The main assumption of this approximation is that 
the external (confining) potential varies smoothly on the 
length scale of the magnetic length, lb — ^/h/ (mujc), 
where m is the effective mass of an electron in a 
GaAs/AlGaAs hetero-structure, and ujc is the cyclotron 
frequency given by tOc = eB/m for the magnetic field 
strength B. At the magnetic field strengths of our in- 
terest, where the average filling factor {D) is around 2, 
i.e. B > ST, /f, is on the order of 10 nanometers, hence 
the TFA is valid. We note that spin degeneracy will not 
be resolved in our calculations. This can be done if the 
cyclotron energy is much larger than the Zeeman energy 
(i.e. effectively we set g = 0). 

In the following, we briefly summarize the self- 
consistent numerical scheme adopted in this work. We 
will assume the 2DES to be located in the plane z = 
85nm with a (surface) number density nci{x, y). We con- 
sider a rectangle of finite extent ax x ay in the aiy-plane, 
with periodic boundary conditions. The (Hartree) con- 
tribution Vff(x,y) to the potential energy of an electron 
caused by the total charge density of the 2DES can be 
written as^ 

2g2 ra^c ray 

VHix,y) = — / dx'dy' K{x,x\y,y') Helix' ,y'), 

1^ Jo JQ 

(2) 

where — e is the electron charge, n an average background 
dielectric constantr^ and the kernel K{x,x' ,y,y') de- 
scribes the solution of Poisson's equation with appropri- 
ate boundary conditions. This kernel can be found in a 
well known text book^. The electron density in turn is 
calculated in the Thomas-Fermi approximation (TFA)^ 

nei(a;,y)= J dE D{E)f{[E + Vix,y) ^ fi*]/kBT), (3) 

with D{E) the relevant (single-particle) density of states 
(DOS), /(s) = [1 + 6"]-^ the Fermi function, and /i"^ the 
electrochemical potential. The total potential energy of 
an electron, V{x,y) = Vcxtix,y) + VH{x,y), differs from 
VH{x,y) by the contribution due to external potentials, 
e.g. the confinement potential generated by the QPC (see 
figure|21), potentials due to the donors etc. The local (but 
nonlinear) TFA is much simpler than the corresponding 
quantum mechanical calculation and yields similar re- 
sults if V{x, y) varies slowly in spaca^, i.e. on a length 
scale much larger than typical quantum lengths such as 
the extent of wave functions or the Fermi wavelength. 

IV. NUMERICAL CALCULATIONS 

The equations lj2J) and (j2J) have to be solved self- 
consistently for a given temperature and magnetic field, 
until convergence is obtained. In our scheme we start 
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with vanishing field and at zero temperature to obtain 
the electrostatic quantities and use these results as an 
initial value for the finite temperature and field calcu- 
lations. For B,T > we start with a relatively high 
temperature and reduce T stepwise in order to obtain a 
good numerical convergence. 



A. Zero magnetic field 

In this section we review the theory of screening in a 
homogeneous 2DES. 

Mesoscopic systems like quantum dots, Hall bars, or 
any edges of quasi-2D electron systems are defined by 
lateral confinement conditions, which lead to an inho- 
mogeneous electron density. An exact treatment of the 
mutual interactions of the electrons in such systems is 
only possible for quantum dots with very few (less than 
ten) electrons. 

The total potential seen by any electron is given 
by the sum of the external potential (describing the 
confinement) and the Hartree potential given by Eq. 
(O, where the electron density in turn is determined 
sclf-consistently by the effective single-particle potential 

Now consider a 2DES in the xy- plane (with vanishing 
thickness) and having the charge density 



n«e*'i""(5(2) (4) 



with q = 27r(nj;/a, nj,/5), where is the q- compo- 
nent of the Fourier transformed electron density. We 
want to obtain the effects of an external perturbation 
5Vext{Y, z), whose Fourier components in the plane z = 
are 5V^^^{Q). This potential induces a charge density (Jn*, 
which in turn leads to an induced potential 



(5) 



that has the tendency to screen the applied external 
potential. Within the TEA, the induced density is re- 
lated to the overall screened potential Vscr by 5nfl = 
-Dl^SV^l^{0), where is the thermodynamic DOS 
defined as Dt ^ J dED{E)^. 

Employing SVscr = SVext + SVind, this yields 



where 



r(q) = l + — 



(6) 



(7) 



is the 2D dielectric function with the Thomas-Fermi mo- 
mentum 



qTF 



(8) 



Then the self-consistent potential at distance \z \ from the 
2DES is 



qTF 

q + qTF 



-e-'l^l<5KL(0), (9) 



i.e., the screening effect of the 2DES decreases exponen- 
tially with \z\. 

In the limit B = 0, T ^ and with Ep = pL*{B ^ 
0,T = 0), Eq. lO reduces to 

n,i{x,y) = Do {Ep - V{x,y)) e{EF - V{x,y)), (10) 

where Dq is the constant DOS for a 2DES given by = 
m/['Kh^). This is a linear relation between V{x,y) and 
nci{x,y) for aU V{x,y) < Ep- 

Now we apply these results to determine the screen- 
ing of a given periodic charge distribution in the plane 
z = 0, which creates an external potential T43;t(r, 0) — 
V'^e^'^'^ in this plane. The self-consistent potential in 
a 2DES then is described by: 



.{r,z)^J2^,Uz)e 



( 



z) = y9e-«^(H- — 



(11) 

The dielectric function e{q) can be expressed in terms 
of the effective Bohr radius = Rfi? / ime^) (for GaAs 
= 9.8 nm), sinceii*^ 2/a^ = 2TTe^Do/R, with q = 
27r/a. We will assume that £{q) ^ 1, so that the TEA is 
valid for i? > 1 T, i.e. < 30 nm. We also note that the 
(7 = component is cancelled by the homogeneous donor 
distribution, assuring overall charge neutrality. 



B. Simple example: Square gate barrier 

We start our discussion by a simple example that 
presents the features of non-linear screening in a 2DES. 
We assume a negatively charged metallic square gate de- 
picted by the white area in the inset of|3i, located at the 
center of a cell that is periodically continued through- 
out the plane (with periods ax = ay = 600nm). The 
square is of size 200nm, and it is kept at the gate poten- 
tial, Vgate = — O.IV. In Eig|21we show the external and 
the screened potential for different separation distances 
of the 2DES and the gate, calculated along the dashed 
line shown in the inset, in the plane of the 2DES. 

In the left panel, the external potential has been plot- 
ted, with the dashed line representing the barrier (gate 
potential) on the surface. We observe that the potential 
profile becomes smooth quickly due to the exponential 
decay of the amplitude of Fourier components at large q 
with increasing a^. 

In contrast, the screened potential displays an interest- 
ing, strong feature close to the edges of the gate [x ~ 200 
and X ^ 400), when the separation distance is relatively 
small (z < 60nm). This is nothing but the manifesta- 
tion of the q dependent screening given in Ea.( lll|l : The 
large q components of the potential remain (almost) unaf- 
fected by screening, whereas the low q (long wavelength) 



5 



0.10- 

> 

0} 

X 0.08- 
N 

o" 

O 

0.06- 



0.02 



-^600nm— - 
















200nm 





I 



^ 0.010 
O 0.005 

o 

CO 

0.000 



b) 





z-Onm 

z=2nm 

— — z=15nm 

z=60nm 
z=85nm 

z=120nm 
z-200nm 




200 400 

X{nm) 



600 



200 400 600 

X(nm) 



FIG. 2: External potential (left panel) seen by a 2DES at dif- 
ferent distances z and the corresponding screened potentials 
(right panel). The separating dielectric material is assumed 
to be GaAs with k = 12.4 and the calculations are done at 
r = OK. 



components are well screened. As a result, wc observe 
sharp peaks near the edge of the gate for small distances 
z, which turn into "shoulders" at larger z. We should 
caution, however, that for z < 60nm the validity of the 
TFA may become questionable, since the potential then 
changes rapidly on the scale of the Fermi wavelength. 

This simple example already demonstrates the strongly 
non-linear behavior of the screening, which can be sum- 
marized as follows: (i) the strongly varying part (high q 
components) of the external potential remains (almost) 
unscreened by the 2DES, but its amplitude decreases 
fast with increasing separation z, whereas (ii) the slowly 
varying part (small q components) is well screened by 
the 2DES, but its amplitude decays much slower for 
large separation distances. Indeed this non-linearity (q- 
dependence of e) leads to pecuhar effects both on electro- 
statics and transport properties of the QPCs, depending 
on the geometry and the structure of the sample. In the 
next section we will look for such effects with regard to 
the QPCs. 



C. Simulation of the QPC 

In this section we will first obtain the bare confinement 
potential created by the QPC for the geometry given in 
figure^ and then go on to discuss the effects of screening. 
The potential generated by such gates can be calculated 
by the scheme proposed by Davies et alm^. 

The model parameters are taken from the relevant ex- 
perimental samples2il4, where the applied gate voltage is 
—0.3V, the width at the tip is about 200nm, and the tip 
separation Ay ^ 300nm. The 2DES is taken to be 85nm 
below the surface. 

We define the QPC using rectangles and polygons 



which are shown in figure ^ as red (dark) and white ar- 
eas. In figure 13 we show the bare confining potential for 
the parameters given above. The electrons are filled up 
to the Fermi energy [Ep ~ 7meV, corresponding to a 
typical electron surface density Uei ~ 1.7 x 10^^ cm^^). 
Using such parameters, the full screening calculation to 
be discussed below will reveal the electrons to be depleted 
beneath the QPC, say at all the dark (blue) regions in 
FigH 

In our numerical simulations, we have mapped the 
unit cell containing the QPC of physical dimensions 
3.3^mxl.8/im to a matrix of 200 by 200 mesh points 
in the absence of a magnetic field and l.l/xmxl.8^m to 
a matrix of 48 by 96 mesh points in the presence, which 
allows us to perform numerical simulations within a rea- 
sonable computation time. With regard to numerical ac- 
curacy, we estimate that, for typical electron densities, 
the mean electron distance, i.e. the Fermi wavelength, 
is larger than 40nm. Hence, the number of mesh points 
considered here allows us to calculate the electron den- 
sity with a good numerical accuracy. We also performed 
calculations for finer meshes and the results do not dif- 
fer quantitatively (at the accuracy of line thicknesses), 
whereas the computational time grows like the square 
of the number of the mesh points. We should also note 
that due to computation time concerns we had to use a 
smaller unit cell in the presence of the magnetic field, 
which yields finite size effects close to the boundaries of 
the sample (e.g. see Fig[7l3). The features observed are, 
in principle, negligible and they tend to disappear when 
the unit cell is taken to be larger and mapped on a larger 
matrix. 

We now discuss the resulting bare and screened poten- 
tial for a realistic QPC defined by surface gates, with a 
tip opening Ay = 300nm. Figure |21 represents the exter- 
nal potential created by the QPC gate structure at the 
surface, calculated in the plane of the 2DES located at 
z = 85nm below the surface, with an applied potential 
—0.3V. In the upper panel we show a 3D plot and a pla- 
nar projection, together with four guide lines, which in- 
dicate the location of the cross-sections through the that 
are displayed in the lower panel. The level of the Fermi 
energy of the system (to be assumed below) is indicated 
in the 3D plot as well. These results have been obtained 
numerically from Eq. The barrier is formed by the 
regions of elevated potential. 

At the first glance one observes that the potential land- 
scape is smoothly varying. This is purely an effect of the 
relatively large distance to the gate, as screening effects 
have not yet been included. For the given Fermi energy 
(obtained from the electron density in the bulk) and the 
tip separation. Ay > 100 nm, the number density of elec- 
trons inside the QPC opening fulfills the validity relation 
of the TFA, i.e. ne;(center)aB ^ 1. At the positions 
where the height of the barrier becomes larger than the 
Fermi energy (light line in the 3D plot and horizontal 
dashed line in the lower panel), the probability to find 
an electron is zero within the TFA. 
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FIG. 3: The bare confinement potential generated by the 
QPC, defined by the polygons shown in figure The color 
scale indicates the strength of the confinement. 



We proceed in our discussion with a comparison of the 
screened potential shown in Fig. 21 to the bare confine- 
ment potential discussed up to now (Fig. 13). The self- 
consistent potential is obtained from the formalism de- 
scribed above for periodic boundary conditions at zero 
temperature and zero magnetic field. The electrons arc 
filled up to the Fermi energy (shown by the gray thick 



line on the surface of the color plot and dashed line in 
the lower panel), such that no electrons can penetrate 
classically into the barrier above those lines. The first 
observation is that the potential profile becomes sharper 
for the screened case and strong variations are observed 
in the vicinity of the QPC. These shouldcr-likc local max- 
ima near the QPC represent the same feature seen in the 
example of the square barrier discussed previously, and 
we have pointed out that they stem from g-dependent, 
non-linear screening. This will become more important 
when we consider a magnetic field, since the local "pin- 
ning" of the Landau levels to the Fermi energy in these 
regions will produce compressible regions surrounded by 
incompressible regions. 

An interesting feature occurs near the opening of the 
QPC, namely a local minimum which is a result of the 
non-linear screening. We point out that somewhat simi- 
lar physics has been found (using spin-dcnsity-functional 
theory^) to lead to the formation of a local bound state 
inside a QPC, which has been related to the "0.7" 
anomaly, linking it with Kondo physics. We believe this 
feature to be a very important result of the self-consistent 
screening calculation, and we note that it may affect 
strongly the transport properties of the QPC both in the 
presence or absence of a magnetic field. We will discuss 
the influence of this local minimum on the formation of 
the incompressible strips in section llV Dl where wc calcu- 
late the density and potential profiles including a strong 
perpendicular magnetic field. 

It is known from the experiments that the interference 
pattern and the transmission properties strongly depend 
on the structure of the QPCs, such as the distance of 
the 2DES from the surface, the applied gate voltage, the 
sharpness and the geometry of the edges, as well as the 
width of the opening of the QPC. The effect of the first 
two parameters can be understood by following the sim- 
ple arguments of linear screening as shown for the square 
gate model: if the distance from the QPC to the 2DES 
increases, the potential profile becomes more and more 
smooth. The screened potential changes linearly with 
the applied gate potential (see Ea. lfTT|) '). The geometric 
parameters have to be adapted to the experiment in ques- 
tion. Note that the shape of the QPCs has already been 
discussed in the literature (see Ref. |2 and references con- 
tained therein). The effect of the size of the QPCs, how- 
ever, has not been considered for large Ay (> lOOnm), 
and we believe this to be an important parameter for the 
interferometer experiments. 

We start our investigation by looking at the opening 
of the QPC with increasing tip separation of the metal 
gates used to define the QPC. In this section, we work 
at zero temperature and magnetic field, with a constant 
bulk electron density. 

In figureElwe depict the self-consistent potential at the 
center of the QPC (j/ = 550nm), while changing the tip 
separation (Ay) between 100 and 500nm. We see that 
for the narrowest separation the potential profile looks 
rather smooth and a minimum is observed at the center. 
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FIG. 5: The screened potential at x = 550nm for five different 
tip distances. Note that the x interval used for the calculation 
smaller than in the previous figure, since we concentrate on 
the bulk structures rather than the edge ones. 



V„(x,y) (eV) 




. 0.006 



J- 2000 \ 
V 3000 



> 



0.010 



— y= Onm 

— y=468nm 
■ ^ y=680nm 

— y=750nm 
y=900nm 
Ep=7.25 meV 



"0.005 



0.000 




1000 1500 2000 2500 3000 

X (nm) 



FIG. 4: The screened potential (upper panel) seen by a 2DES 
at 85nm below the surface and some characteristic cuts along 
the X- axis, together with an indication of the Fermi level Ef 
(lower panel). The color scale represents the strength of the 
potential, and the cross-sections are indicated by the same 
line code as in Fig. |3 



If we increase Ay(^ 300nm) we see that the screening be- 
comes stronger, leading to more pronounced shoulders on 
the sides and a deeper minimum at the center. For even 
larger separations (Ay > 300nm) a local maximum starts 
to develop at the center, since the electrostatic potential 
energy is no longer strong enough to repel the electrons 
from this region. Basically all the non-linear features 
observed result from the competition between the gate 
potential, which simply repels the electrons, and the mu- 
tual Coulomb interaction, i.e. the Hartree potential. It 
is obvious that for narrower tip separations only a few 
electrons will remain inside the QPC opening and there- 
fore TFA type approximations will not be justified any 
longer. 

Summarizing this section, we have determined the 
screened potential profile in a realistic QPC geometry, 
pointing out features resulting from non-linear screen- 
ing. We have observed that a local extremum occurs 
at the center of the QPC, and have traced the depen- 
dence on the width Ay between the QPC tips. These 
features, as mentioned before, become more interesting 
if a magnetic field is also taken into account, where they 
lead to stronger spatial inhomogeneities in the electron 
distribution. Our next step is thus to include a strong 
quantizing perpendicular magnetic field and examine the 
distribution of the incompressible strips where the im- 
posed external current is confined^SiS^. 



D. Finite temperature and Magnetic field 



Once the initial values of the screened potential and the 
electron distribution have been obtained for T = 0, B = 
0, using the scheme described above, one can calculate 
these quantities for finite field and temperature as fol- 
lows: replace the zero temperature Fermi function with 
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FIG. 6: [a]-[d] Color-coded plot of the local filling factor ver- 
sus position {x, y) for a square sample of width = ay = 
0.8 /xm; white indicates ^{x, y) — 2. The average density is 
taken to be 3.0 ■ lO" cm"^; kBT/Ep = 0.02. [e] A sketch of 
the Hall resistance as a function of magnetic field 



the finite temperature one and insert the bare Landau 
DOS 



Db{E) = —Y^5{E~ En), En = hLo,{n+l/2) (12) 

n=0 

into equation Q instead of Dq. In our numerical scheme 
we first start with relatively high temperatures (i.e. a 
smooth Fermi function) and then decrease the tempera- 
ture slowly until the desired temperature is reached. A 
Newton- Raphson method is used for the iteration process 
and at every iteration step the electro-chemical potential 
is checked to be constant. 

Before proceeding with the investigation of the QPC 
geometry at B > 0, we would like to make clear the re- 
lation between the quantum Hall plateaus and the exis- 
tence of the incompressible strips following the arguments 
of Siddiki and GerhardtsS^. Fig. presents the local fill- 
ing factors of a relatively small Hall bar. together with an 
illustrative Hall resistance curve. At the high magnetic 
field side ((b), iy{0,0) < 2) there are no incompressible 



strips, thus the system is out of the Hall plateau. When 
approaching from the high B side to the plateau a single 
incompressible strip at the center develops. When the 
width of this strip becomes larger than the Fermi wave- 
length, the system is in the quantum Hall state ((d), 
1/(0, 0) = 2). If we decrease the field strength further the 
center incompressible strip splits into two and moves to- 
ward the edges ((c) !^(0,0) > 2). As long as the widths 
of these strips are larger or comparable with the Fermi 
wavelength the system remains in the plateau. This is 
the regime in which an interferometer may be realized. 
Further decreasing the magnetic field leads to narrower 
incompressible strips which finally disappear if the widths 
of them become smaller than the average electron dis- 
tance. Then the system leaves the quantized plateau. 
The distribution of the incompressible strips and the on- 
set of the plateaus, of course, depends on the disorder 
potentials^ and the physical size of the sample. However, 
the experiments considered here are done using narrow 
and high mobility structures, thus the above scheme will 
cover the experimental parameters. 

In this subsection we present some of our results ob- 
tained within the TFA using periodic boundary condi- 
tions, considering two different tip separations, while 
sweeping the magnetic field. First we will fix the gate 
potential to Vqpc = —0.3V and sweep the magnetic field 
for Ay = lOOnm, while keeping the electron number den- 
sity, i.e. the Fermi energy, constant. Second we examine 
the potential profile for Ay = 300nm and comment on 
the possible effects on the coherent transport properties. 

In figure[71wc plot the local filling factor (i.e. the nor- 
malized density) distribution of the 2DES projected on 
the ccy-plane, together with the same quantity for some 
selected values of y, at average filling factor {D) one. 
From the y = Onm curve (solid lines) in FigEIb], one 
can see that the electrons beneath the QPC arc depleted 
(shaded, dark (blue) regions) (300 < x < 800nm), while 
the electron density reaches finite values while approach- 
ing the opening of the QPC {y ~ 850nm). At D ~ 1 
one does not observe any incompressible regions, since 
the Fermi energy is pinned to the lowest Landau level. 
Hence the electron distribution is rather smooth and the 
current distribution will just be proportional to the num- 
ber of electrons, similar to the Drude approach. For this 
case the external potential is screened almost perfectly 
and the self-consistent potential is almost flat, thus one 
can assume that the corresponding local wave functions 
are very similar to the ground state Landau wave func- 
tions. 

The first incompressible region occurs when the Fermi 
energy falls in the gap between two low-lying Landau 
levels. Then the electrons exhibit a constant density 
and thus cannot screen the external potential. In fig- 
ure ©a, we show the electron distribution for D = 1.1. 
The black regions denote a local density corresponding 
to filling factor i/ = 2, which does not percolate from 
the left side of the sample (which we might identify with 
the source) to the right side (drain). Here one can see 
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FIG. 7: [a] The top view of the local filling factor, v{x,y), 
distribution of the 2DES, for average filling factor one, in the 
plane located at 2 = 85nm below the surface, at the "de- 
fault" temperature, kT/huJc = 1/50, which will be used in 
all subsequent plot. The color scale depicts the density of 
electrons, whereas the dark shaded areas indicate the elec- 
tron depleted regions, [b] Side view of the local filling factor 
for y = Onm(solid line), 468nm (dashed line), 660nm (dot- 
ted line), 750nm (dash-dotted line) and 900nm (dash-dotted- 
dashed line) . The horizontal lines in [a] shows the positions 
of the cuts in [b], with the same line code. Note that the den- 
sity has local minima at large and small x, which are finite 
size efi'ects mentioned in the text. 



well developed incompressible puddles, at the regions 
150nm< X < 250nm, Onm< y < 450nm (and four other 
symmetric ones) and two smaller puddles at the entrance 
of the QPC. These structures will remain unchanged even 
if one considers a larger unit cell, since they manifest the 
q- dependency, i.e. the rapid oscillations of the Fourier 
transform of the confining potential of the QPC. 

In these regions the self-consistent potential exhibits 



a finite slope. Accordingly the wave functions will be 
shifted and squeezed, i.e. they are now superpositions 
of a few high order Landau wave functions with re- 
normalized center coordinates. This behavior has been 
showni2Si2& for the translationally invariant model. Here 
we did not include the finite extent of the wave functions, 
here to avoid lengthy numerical calculations. 

The incompressible regions shift their positions on the 
xy plane depending on the strength and the profile of 
the confining potential. In figure JHJb we show the filling 
factor distribution where the bulk filling factor is almost 
two. We see that four incompressible strips arc formed 
near the QPC. However the QPC opening remains in 
a compressible state, with local filling factor less than 
two, where we expect that the self-consistent potential 
is essentially flat. Further increasing the average filling 
factor, we observe that the bulk becomes completely com- 
pressible and two incompressible strips are formed near 
the QPC which percolate from bottom to top, creating a 
potential barrier with a height of hujc, see Fig. (jS);). For 
even higher filling factors, they merge at the center of 
the QPC (Fig.|Hli). In that case, the potential within the 
QPC will then no longer be flat, due to poor screening. 
We should also note that for a small width Ay of the 
QPC opening, merging of the incompressible strips will 
happen only in a very narrow B interval, and a quan- 
titative evaluation within our TFA can not be always 
satisfactory, as the number of electrons inside the QPC 
becomes too low. Further decreasing the field strength 
(increasing the average filling factor) results in two sep- 
arate incompressible strips winding around the opposite 
gates making up the QPC, as shown in Fig. Thus, 
dissipationless transport through the QPC, with a quan- 
tized conductance, becomes possible. At the lowest field 
values considered in this figure, we see that the innermost 
incompressible strips (with 1/ — 2) become smaller than 
the Fermi wavelength and thus they essentially disappear 
and no longer affect the transport properties. This point 
has been discussed in detail in a recent work by Siddiki 
et alm^. The scheme described above now starts to re- 
peat, but with incompressible strips having a local filling 
factor of 4. 

We now discuss the effects of increasing the separation 
parameter, which we choose to be Ay = 300nm in figure 
© . At the strongest magnetic field ^j^) , only very small 
regions are incompressible and the electron distribution is 
similar to Fig|SJi, where the incompressible regions result 
from local unpinning of the Fermi energy from the low- 
est Landau level due to dependent screening, i.c the 
shoulder-like variation of the potential near the QPC dis- 
cussed earlier. By decreasing B, an interesting structure 
is observed at the center of the QPC: an incompress- 
ible island. In figure ^jp, we have tuned the magnetic 
field such that the bulk of the 2DES is incompressible, 
meanwhile the entrance to the QPC remains compress- 
ible. The strong variation of the self-consistent potential 
at the center of the QPC can generate a pronounced effect 
on the current passing through the QPC (see figure [TUI 
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FIG. 8: The local filling factor distribution for different av- 
erage filling factors {v), which is defined by the number of 
the electrons in the unit cell, [a] i> = 1.1 [b] 1.2 [c] 1.4 [d] 
1.6 [e] 1.8 [f] 2.2. The color scale depicts the local electron 
concentration, whereas the abrupt colors indicate the even- 
integer filling factors, i.e. incompressible strips, (black for 
v{x,y) = 2, magenta for i/{x,y) = 4. 



and the related text). For a lower magnetic field strength 
the distribution of the incompressible region is just the 
opposite (c). Now we see a large compressible puddle at 
the center, surrounded by incompressible regions, which 
can percolate from source to drain. Coherent, dissipa- 
tionless transport can be expected in this case. Further 
decreasing the magnetic field we observe that the struc- 
ture is smeared out and the tip region becomes com- 
pressible, nevertheless there are two large incompressible 
regions close to the entrance of the opening. The two 
incompressible strips wind around the gates, as shown 
in Fig. Ol. Finally, a scheme similar to that observed 
earlier in figure |Hli-e is also seen now, for Ay = 300nm. 

Another remark which we would like to make concerns 
the edge profile of the sample itself and of the QPC. It has 
been shown both cxperimentallji^ and thcorcticalljiSSi^ 
that for an (almost) infinite potential barrier at the edges 
of the sample, the Chklovskii edge state picture 
breaks down, i.e. no incompressible strips near the edge 
can be observed. Meanwhile for smoothly varying edge 
potential profiles many incompressible strips are present, 
if the bulk filling factor is larger than two (for spinless 
electrons: four). We believe that, within the MZI setup 
both of these edge potential profiles might co-exist. At 
the edge regions of the sample, where lateral confine- 
ment is defined by physical etching, the potential profile 
differs from of the one generated by the top gates, due 
to different separation thicknesses and also lateral sur- 
face charges generated by etching. In principle gate and 
etching defined edges impose different boundary condi- 
tions, and the effects on screening at a 2DES have been 
discussed before24. These two profiles will certainly af- 



fect the group velocity, since the slope of the potential 
depends on the (lateral) boundary conditions. Follow- 
ing the arguments of Ref . [2^l25| . which essentially pre- 
dict that the dissipative current is confined to the in- 
compressible strips, the widths of these strips will also 
define the slope, hence the velocity of the electrons will 
be determined by the edge profile. The velocity of the 
edge electrons were investigated experimentallji^ and the 
magnetic field dependency was reported as B^^^. There 
it was noted that a self-consistent treatment is necessary 
to understand their findings, which we would like to dis- 
cuss in a future publication. 

The important features to note in these results are (i) 
in general, electron-electron interactions have a remark- 
able effect, leading to the formation of a local extremum 
in the potential at the center of the QPC, which even at 
low electron densities seems to be well described by the 
TFA; (ii) the narrow compressible/incompressible strips 
formed near the QPC are a direct consequence of the 

dependent screening. 



E. Comments on coherent transport 

A complete calculation of coherent transport requires 
a deeper analysis of the wave functions and is beyond 
the scope of this work, which has been devoted to self- 
consistent realistic calculations of the potential and den- 
sity profiles. In principle, one can follow the arguments of 
the well developed recursive Green's function technique^ 
in the absence of magnetic field and the method devel- 
oped recently even in the presence of a strong field^. 

Instead we would like to examine the potential dis- 
tribution across the QPC and comment on the possible 
effects of interaction on the wave functions, and thereby 
(indirectly) on transport. In figure [TUI we depict the po- 
tential profile across the QPC for the parameters used to 
obtain figure |S1 As expected for D = 1.0 (dashed (red) 
line) the 2DES is "quasi" metallic, hence the external 
potential is perfectly screened, and the wave functions 
are left almost unchanged. The two incompressible is- 
lands seen at the entrance of the QPC in Fig|Hli lead to a 
minor variation of the screened potential at a; = 300nm 
and X = 800nm, depicted by the solid (black) line for 
D ~ l.l. A drastic change is observed when the bulk be- 
comes incompressible {D = 1.2) and the opening remains 
compressible: Now the 2DES cannot screen the exter- 
nal potential near the openings of the QPC, where we 
see a strong variation. The strong perpendicular mag- 
netic field changes the potential profile near the QPC 
via forming incompressible strips, and local minima are 
observed at the entrance and the exit. In these regions 
the electrons are strongly localized and the wave func- 
tions are squeezed. The situation is rather the opposite 
for D ~ 1.4, where two incompressible strips located near 
the QPC, formed due to g-depcndent screening, merge 
at the opening. One observes a barrier with the height 
of fkUc, which essentially is a direct consequence of the 
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FIG. 9: The local filling factor distribution for different av- 
erage filling factors (i?), for a tip separation 6y = 300 nm. 
Note that, the number of electrons in the unit cell is changed, 
since the depleted areas are larger than of the previous fig- 
ure, [a] v = 1.14 [b] 1.2 [c] 1.34 [d] 1.4, [e] 1.6 and [g] 
3.1. The color scale depicts the local electron concentration, 
whereas the high-contrast color regions indicate the even- 
integer filling factors, i.e. incompressible strips, (black for 
v{x,y) — 2, magenta for v{x,y) = 4). The calculations 
are done at ksT/huJc = 1/50 for an average electron density 
1.7 X 10""cm-^ 



incompressible strip at the center and electrons have to 
overcome this barrier. Further decreasing the magnetic 
field smears out the barrier gradually, until the system 
becomes completely compressible and wc arc back in the 
case of figurelTUh (also with regard to the transport prop- 
erties). 



V. SUMMARY 

The study that motivated the present authors 
was the Quantum Hall effect based Mach-Zender 
interferometeB^ii. There are puzzles in the experiment 
for which it is not obvious (at least to us and some oth- 
ers) how they could be explained using scattering the- 
ory. Therefore we probably need to take into account 
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FIG. 10: The self-consistent potential profile across the QPC, 
plotted for characteristic values of the average filling factor. 
Calculations are done at the default temperature and electron 
density. 



interactions more seriously, including correlation effects 
in non-equilibrium transport that will not be part of scat- 
tering theory. A first step towards that goal then is the 
self-consistent mean-field equilibrium calculation which 
we have done. 

At the moment we do not have to offer a non- 
equilibrium transport theory based on these equilibrium 
calculations. These effects may include decoherence 
due to potential fluctuations brought about by electron- 
electron or electron-phonon interactions (together with 
other noise sources). A more detailed understanding of 
electron-electron interactions in this setup, as well as of 
those features of the interferometer that are specific to 
the physics of the Quantum Hall effect, hinges on an anal- 
ysis of the self-consistent static potential landscape near 
the QPCs, which represent the most crucial components 
of the setup. 

Therefore, in this work, we have taken into account 
the electron-electron interaction within the TFA, consid- 
ering realistic geometries of QPCs, calculating the self- 
consistent potential and electronic density profiles, and 
commenting on possible effects on transport. 

The outcome of our model calculations can be sum- 
marized as follows: (i) We have obtained the electro- 
static potential profile for the QPC geometries used in 
the experiments by solving the Laplace equation semi- 
analytically. (ii) We have demonstrated for a simple 
square well barrier that the screened potential in a 2DES, 
even in the absence of a magnetic field, is strongly depen- 
dent on the initial potential profile and on the distance 
between gates and 2DES. (iii) The screened potential has 
been calculated within the TFA for a QPC at vanishing 
field, where we have observed two interesting features: 
a local extremum at the center of the QPC and strong 
shoulder-like variations near the QPC. (iv) In the pres- 
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ence of a magnetic field, the formation and the evolu- 
tion of the incompressible regions has been examined and 
three different cases have been observed: (a) the system 
is completely compressible, (b) An incompressible re- 
gion and/or strip, which does not percolate from source 
to drain, generates a local extremum at the entrance/exit 
of the QPC (c) The center of the QPC becomes incom- 
pressible, with or without a compressible island, hence 
the incompressible strip percolates from source to drain. 
We note that the local minimum found at the center of 



the QPC for certain tip separations, being a clear inter- 
action effect, coincides with the findings of Hirose et al 
0. 
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In this work we implement the self-consistent Thomas-Fermi-Poisson approach to a homogeneous 
two dimensional electron system (2DES). We compute the electrostatic potential produced inside a 
semiconductor structure by a quantum-point-contact (QPC) placed at the surface of the semicon- 
ductor and biased with appropriate voltages. The model is based on a semi-analytical solution of 
the Laplace equation. Starting from the calculated confining potential, the self-consistent (screened) 
potential and the electron densities are calculated for finite temperature and magnetic field. We ob- 
serve that there are mainly three characteristic rearrangements of the incompressible "edge" states, 
which will determine the current distribution near a QPC. 

PACS numbers: 73.20.-r, 73.50.Jt, 71.70.Di 



I. 



INTRODUCTION 



A quantum point contact (QPC) is constructed by 
geometric or electrostatic confinement of a two dimen- 
sional electron system (2DES). The conductance through 
them is quantized^^ and they play a crucial role in the 
field of mesoscopic quantum transport. Their properties 
have been investigated in a wide variety of experiments, 
which include the observation of the 0.7 anomalji^, 
quantum dots coupled to QPCs^, Quantum-Hall effect 
(QHE) based Mach-ZenderiSii and Aharonov-Bohm in- 
terferometers. This has lead to extensive investigations 
of the electrostatic and transport properties of QPCs, 
both with and without a quantizing magnetic field. Many 
different techniques have been used to find the elec- 
tronic density distribution near a QPC, ranging from nu- 
merical Poisson-Schrodinger solutions^ to spin-density- 
functional theorji^ and phenomenological approachesiS. 
It has been possible to treat realistic samples mostly 
only within simplified electrostatic calculations, neglect- 
ing screening effects. On the other hand, when including 
interactions the calculations become more complicated, 
thus one usually sacrifices handling realistic geometries. 

Recent experiments have succeeded in developing and 
analyzing a QHE based electronic Mach-Zcnder inter- 
ferometer (MZI)^, making use of the integer QHE edge 
stateai as single-channel chiral quantum wires. A key 
element of these experiments are the QPCs, which play 
the role of the beam splitters of the optical setup. The 
electrostatic potential and electronic density distribu- 
tions in and near the QPCs play an important role in 
understanding the rearrangement of the edge states in- 
volved. Moreover, the electron-electron interaction has 
been proposed^ as one of the origins of dephasing in such 
an electronic MZI, such that a self-consistent calculation 
of the electrostatic potential may also be viewed as a first 
step towards a quantitative understanding of this issue. 
So far, the theoretical description of dephasing in the 
electronic MZI via classicaliiiiSiiii^ or quantum noise 
fieldsi^ii^ and other approachesii has focused on features 



supposed to be independent of its specific realization (see 
Ref. ^3 foi' * recent review). However, a more detailed 
analysis of the QHE related physics, taking account of 
interaction effects, will certainly be needed for a direct 
comparison with experimental data. In this paper, we 
will provide a detailed numerical analysis of the electro- 
statics of QPCs in the integer QHE, assuming geometries 
adapted to those used in the MZI experiment. Our work 
will produce the electron density and electrostatic poten- 
tial, based on the self-consistent Thomas-Fermi-Poisson 
approximation, to which we refer as TEA in the following. 

We would like to point out the following observation 
regarding the Mach-Zchnder experiment, where a yet- 
unexplained beating pattern observed in the visibility 
(interference contrast) as a function of bias voltage was 
surprisingly found to have a period independent of the 
length of the interferometer arms. Such a result would 
seem less surprising if all the relevant interaction physics 
leading to the beating pattern were actually taking place 
in the vicinity of the QPC. This provides strong encour- 
agement for future more detailed work on the coherent 
transport properties of these QPCs. 

Although it has been more than two decades since the 
discovery of the quantized Hall effectiS, the microscopic 
picture of current distribution in the sample and the in- 
terplay of the current distribution with the Hall plateaus 
is still under debate. In recent experiments, the Hall 
potential distribution and the local electronic compress- 
ibility have been investigated in a Hall bar geometry by 
a low-temperature scanning force microscopeSS and by 
a single-electron-transistoriSi, respectively. This has mo- 
tivated theoretical^^ work, where a self-consistent TEA 
calculation has been used to obtain electrostatic quanti- 
ties. 

Self-consistent screening calculations show that the 
2DES contains two different kinds of regions, namely the 
quasi-metallic compressible and quasi-insulating incom- 
pressible regions22i24. The electron distribution within 
the Hall bar depends on the "pinning" of the Fermi level 
to highly degenerate Landau levels. Wherever the Fermi 
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level lies within a Landau level with its high density of 
states (DOS), the system is known to be compressible 
(leading to screening and correspondingly to a flat po- 
tential profile), otherwise it is incompressible, with a con- 
stant electron density and, in general, a spatially vary- 
ing potential due to the absence of screening. Moreover, 
based on these results for the potential and density dis- 
tributions, one may employ a local version of Ohm's law 
(together with Maxwell's equations and an appropriate 
model for the conductivity tensor) to calculate the cur- 
rent distribution, imposing a given overall external cur- 
rent for the in-plane geometry . These results are mostly 
consistent with experiments except that within the self- 
consistent TFA one obtains an incompressible strip (IS) 
for a large interval of magnetic field values which leads 
to coexistence of several IS's with different local filling 
factors. Recently, this theory has been improved in two 
aspects2S*2&: (i) the finite extent of the wave functions 
was taken into account in obtaining electrostatic quan- 
tities (rather than using delta functions), (ii) the find- 
ings of the full Hartree calculations were simulated by 
a simple averaging of the local conductivities over the 
Fermi wave length , thereby relaxing the strict locality of 
Ohm's law for realistic sample sizes. A very important 
outcome of this model is that there can exist only one 
incompressible edge state at one side of the sample for a 
given magnetic field value. Indeed this is differing dras- 
tically from the Chklovskii-Shklovskii-Glazman (CSG)^ 
and the Landauer-Biittikei— picture, where more than 
one edge state can exist and is necessary to "explain" 
the QHE. In the CSG scheme a non-self-consistent TFA 
(which is called the "electrostatic approximation") was 
used. However, it is clear that if the widths of the IS's 
(where the potential variation is observed) become com- 
parable with the magnetic length, the TFA is not valid, 
thus the results obtained within this model are not reli- 
able any more. In principle, similar results to Ref. p5l | 
were reported by T. Suzuki and T. Ando2^, quite some 
time ago and recently by S. Ihnatsenka and I. V. Zo- 
zoulenko^ in the context of spin-density-functional the- 
ory. With the improvements on the self-consistent TFA 
mentioned above, together with taking into account the 
disorder potential^ and using the self-consistent Born 
approximation^ to calculate the local conductivity ten- 
sor, one obtains well developed Hall plateaus, with the 
longitudinal resistivity vanishing to a very high accu- 
racy, and one is also able to represent correctly the inter- 
plateau transition regions. Wherever one observes an IS, 
the longitudinal conductivity becomes zero, and as a con- 
sequence also the corresponding local resistance (and the 
total resistance) vanishes. Thus, according to Ohm's law, 
the current flows through the incompressible region. In 
addition, the Hall conductance becomes equal to the local 
value of the quantized conductance. Finally all the three 
experimentally observed^ qualitatively different regimes 
of how the Hall potential drops across the sample have 
been reproduced theoretically without artifacts of the 
TFA2S. The model described above has also been suc- 



cessfully applied to an electron-electron bilayer system^^ 
and provided a qualitative explanation^i of the magneto- 
resistance hysteresis that has been reported recentlji^Si^. 
For all of these reasons, we feel confident in applying this 
theory to our analysis of the MZI setup. 

Motivated by the experimental and theoretical findings 
ascertaining the importance of the interaction effects in 
the integer Quantum Hall regime, in this work we will 
show that the mutual Coulomb interaction between the 
electrons leads to interesting non-linear phenomena in 
the potential and electron distribution in close proximity 
of the QPCs. Based on the self-consistent Thomas-Fermi- 
Poisson approximation (TFA), we will consider realistic 
QPC geometries and examine the distribution of the in- 
compressible regions depending on the field strength and 
sample parameters. 

The rest of this paper is organized as follows: In SecHII 
the electrostatic potential produced by an arbitrary sur- 
face gate will be discussed, by solving the Laplace equa- 
tion without screening effects. In Sec. IIIII we review the 
TFA in a 2DES. In Sec.|lVlwe will first present the well 
known general results of the TFA for a homogeneous 
2DES at zero magnetic field B and zero temperature, 
and we will investigate the electron density and electro- 
static potential profiles of a (i) simple square gate geom- 
etry and (ii) a generic QPC, before (iii) systematically 
investigating the positions of the incompressible strips 
depending on magnetic field and geometric parameters. 
We conclude with a discussion in Sec. 



II. ELECTROSTATICS OF THE GATES 

As mentioned in the introduction, there is a tradeoff 
between simulating realistic QPC geometries and includ- 
ing the interaction effects within a reasonable approxima- 
tion. In this paper, we present an intermediate approach, 
which considers realistic QPC structures but interactions 
of the electrons are handled within a Thomas-Fermi ap- 
proximation (TFA), which is valid for relatively "large" 
QPCs (> 100 nm). One can obtain, in a semi-analytical 
fashion, the electrostatic potential generated by an arbi- 
trary metallic gate at the surface by solving the Laplace 
equation for the given boundary conditions. Afterwards, 
it is possible to obtain the electron and potential distri- 
bution in the 2DES, within the TFA, both for vanishing 
and finite magnetic fields (B), and at low temperatures 
at B > 0. 

Here we briefly summarize the semi-analytical model 
developed by J.Davies and co-workers^. The aim of this 
section is to calculate the electrostatic potential on a 
plane at some position z below the surface of the semi- 
conductor, which is partially covered by a patterned gate. 
The surface occupies the 2 = plane and z is measured 
into the material. The un-patterned surface is taken to 
be pinned so we can set the potential Vup{r, 0) = there, 
with Vgate (r,0) = Vg on the gate. We use lower-case 
letters like r = (x, y) to denote two-dimensional vectors 
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FIG. 1: The image of the QPC (gray scale). The polygons 
are used to define the gates on the two-dimensional mesh. 
The 2DES resides under the dark (grey) regions, with a bulk 
electron density of 1.7 x 10~^^cm~^ (see Ref. due to mem- 
ory limitations the quality of the figure is reduced. The white 
line scales 200nm). 



with the corresponding upper-case letters for three di- 
mensional vectors like R = {x,y,z) = (r, z). Thus the 
problem is to find a solution, VextC^), to the Laplace 
equation V^Veit = 0, given the value on the plane 
z ^ 0, and subject to the further boundary condition 
dVext/dz — > as z — > oo. One route is to start by mak- 
ing a two-dimensional Fourier transform from Vexti^i^) 
to Vext (q, 0) . The dependence on z is a decaying exponen- 
tial to satisfy Laplaces equation and the boundary con- 
dition at z = oo: Vext{(\, z) = Vext{(\, 0) exp(-|qz|). This 
multiplication of the Fourier transform is equivalent to a 
convolution in real space. Taking the two-dimensional in- 
verse Fourier transform of cxp(— |gz|) leads to the general 
result. 

14xt(r,z) = i / l"^' V.ateiy' ,Q)dv' , 

Kj 27r(z^ + |r — r'l"')-^/"' 

(1) 

where k is the dielectric constant of the considered 
hctcro-structure. Now one can evaluate the potential 
in the plane of the 2DES, z = d, for a given gate and 
potential distribution on the surface. The derivation of 
some important shapes like triangle, rectangle and poly- 
gons is provided in the work cited above, which has been 
successfully applied to quantum dot systems^. For our 
geometry, we will use the result for the polygons. 



III. ELECTRON-ELECTRON INTERACTION: 
THOMAS-FERMI-POISSON APPROXIMATION 

The main assumption of this approximation is that 
the external (confining) potential varies smoothly on the 
length scale of the magnetic length, lb — ^/h/ (mujc), 
where m is the effective mass of an electron in a 
GaAs/AlGaAs hetero-structure, and ujc is the cyclotron 
frequency given by tOc = eB/m for the magnetic field 
strength B. At the magnetic field strengths of our in- 
terest, where the average filling factor {D) is around 2, 
i.e. B > ST, /f, is on the order of 10 nanometers, hence 
the TFA is valid. We note that spin degeneracy will not 
be resolved in our calculations. This can be done if the 
cyclotron energy is much larger than the Zeeman energy 
(i.e. effectively we set g = 0). 

In the following, we briefly summarize the self- 
consistent numerical scheme adopted in this work. We 
will assume the 2DES to be located in the plane z = 
85nm with a (surface) number density nci{x, y). We con- 
sider a rectangle of finite extent ax x ay in the aiy-plane, 
with periodic boundary conditions. The (Hartree) con- 
tribution Vff(x,y) to the potential energy of an electron 
caused by the total charge density of the 2DES can be 
written as^ 

2g2 ra^c ray 

VHix,y) = — / dx'dy' K{x,x\y,y') Helix' ,y'), 

1^ Jo JQ 

(2) 

where — e is the electron charge, n an average background 
dielectric constantr^ and the kernel K{x,x' ,y,y') de- 
scribes the solution of Poisson's equation with appropri- 
ate boundary conditions. This kernel can be found in a 
well known text book^. The electron density in turn is 
calculated in the Thomas-Fermi approximation (TFA)^ 

nei(a;,y)= J dE D{E)f{[E + Vix,y) ^ fi*]/kBT), (3) 

with D{E) the relevant (single-particle) density of states 
(DOS), /(s) = [1 + 6"]-^ the Fermi function, and /i"^ the 
electrochemical potential. The total potential energy of 
an electron, V{x,y) = Vcxtix,y) + VH{x,y), differs from 
VH{x,y) by the contribution due to external potentials, 
e.g. the confinement potential generated by the QPC (see 
figure|21), potentials due to the donors etc. The local (but 
nonlinear) TFA is much simpler than the corresponding 
quantum mechanical calculation and yields similar re- 
sults if V{x, y) varies slowly in spaca^, i.e. on a length 
scale much larger than typical quantum lengths such as 
the extent of wave functions or the Fermi wavelength. 

IV. NUMERICAL CALCULATIONS 

The equations lj2J) and (j2J) have to be solved self- 
consistently for a given temperature and magnetic field, 
until convergence is obtained. In our scheme we start 
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with vanishing field and at zero temperature to obtain 
the electrostatic quantities and use these results as an 
initial value for the finite temperature and field calcu- 
lations. For B,T > we start with a relatively high 
temperature and reduce T stepwise in order to obtain a 
good numerical convergence. 



A. Zero magnetic field 

In this section we review the theory of screening in a 
homogeneous 2DES. 

Mesoscopic systems like quantum dots, Hall bars, or 
any edges of quasi-2D electron systems are defined by 
lateral confinement conditions, which lead to an inho- 
mogeneous electron density. An exact treatment of the 
mutual interactions of the electrons in such systems is 
only possible for quantum dots with very few (less than 
ten) electrons. 

The total potential seen by any electron is given 
by the sum of the external potential (describing the 
confinement) and the Hartree potential given by Eq. 
(O, where the electron density in turn is determined 
sclf-consistently by the effective single-particle potential 

Now consider a 2DES in the xy- plane (with vanishing 
thickness) and having the charge density 



n«e*'i""(5(2) (4) 



with q = 27r(nj;/a, nj,/5), where is the q- compo- 
nent of the Fourier transformed electron density. We 
want to obtain the effects of an external perturbation 
5Vext{Y, z), whose Fourier components in the plane z = 
are 5V^^^{Q). This potential induces a charge density (Jn*, 
which in turn leads to an induced potential 



(5) 



that has the tendency to screen the applied external 
potential. Within the TEA, the induced density is re- 
lated to the overall screened potential Vscr by 5nfl = 
-Dl^SV^l^{0), where is the thermodynamic DOS 
defined as Dt ^ J dED{E)^. 

Employing SVscr = SVext + SVind, this yields 



where 



r(q) = l + — 



(6) 



(7) 



is the 2D dielectric function with the Thomas-Fermi mo- 
mentum 



qTF 



(8) 



Then the self-consistent potential at distance \z \ from the 
2DES is 



qTF 

q + qTF 



-e-'l^l<5KL(0), (9) 



i.e., the screening effect of the 2DES decreases exponen- 
tially with \z\. 

In the limit B = 0, T ^ and with Ep = pL*{B ^ 
0,T = 0), Eq. lO reduces to 

n,i{x,y) = Do {Ep - V{x,y)) e{EF - V{x,y)), (10) 

where Dq is the constant DOS for a 2DES given by = 
m/['Kh^). This is a linear relation between V{x,y) and 
nci{x,y) for aU V{x,y) < Ep- 

Now we apply these results to determine the screen- 
ing of a given periodic charge distribution in the plane 
z = 0, which creates an external potential T43;t(r, 0) — 
V'^e^'^'^ in this plane. The self-consistent potential in 
a 2DES then is described by: 



.{r,z)^J2^,Uz)e 



( 



z) = y9e-«^(H- — 



(11) 

The dielectric function e{q) can be expressed in terms 
of the effective Bohr radius = Rfi? / ime^) (for GaAs 
= 9.8 nm), sinceii*^ 2/a^ = 2TTe^Do/R, with q = 
27r/a. We will assume that £{q) ^ 1, so that the TEA is 
valid for i? > 1 T, i.e. < 30 nm. We also note that the 
(7 = component is cancelled by the homogeneous donor 
distribution, assuring overall charge neutrality. 



B. Simple example: Square gate barrier 

We start our discussion by a simple example that 
presents the features of non-linear screening in a 2DES. 
We assume a negatively charged metallic square gate de- 
picted by the white area in the inset of|3i, located at the 
center of a cell that is periodically continued through- 
out the plane (with periods ax = ay = 600nm). The 
square is of size 200nm, and it is kept at the gate poten- 
tial, Vgate = — O.IV. In Eig|21we show the external and 
the screened potential for different separation distances 
of the 2DES and the gate, calculated along the dashed 
line shown in the inset, in the plane of the 2DES. 

In the left panel, the external potential has been plot- 
ted, with the dashed line representing the barrier (gate 
potential) on the surface. We observe that the potential 
profile becomes smooth quickly due to the exponential 
decay of the amplitude of Fourier components at large q 
with increasing a^. 

In contrast, the screened potential displays an interest- 
ing, strong feature close to the edges of the gate [x ~ 200 
and X ^ 400), when the separation distance is relatively 
small (z < 60nm). This is nothing but the manifesta- 
tion of the q dependent screening given in Ea.( lll|l : The 
large q components of the potential remain (almost) unaf- 
fected by screening, whereas the low q (long wavelength) 
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FIG. 2: External potential (left panel) seen by a 2DES at dif- 
ferent distances z and the corresponding screened potentials 
(right panel). The separating dielectric material is assumed 
to be GaAs with k = 12.4 and the calculations are done at 
r = OK. 



components are well screened. As a result, wc observe 
sharp peaks near the edge of the gate for small distances 
z, which turn into "shoulders" at larger z. We should 
caution, however, that for z < 60nm the validity of the 
TFA may become questionable, since the potential then 
changes rapidly on the scale of the Fermi wavelength. 

This simple example already demonstrates the strongly 
non-linear behavior of the screening, which can be sum- 
marized as follows: (i) the strongly varying part (high q 
components) of the external potential remains (almost) 
unscreened by the 2DES, but its amplitude decreases 
fast with increasing separation z, whereas (ii) the slowly 
varying part (small q components) is well screened by 
the 2DES, but its amplitude decays much slower for 
large separation distances. Indeed this non-linearity (q- 
dependence of e) leads to pecuhar effects both on electro- 
statics and transport properties of the QPCs, depending 
on the geometry and the structure of the sample. In the 
next section we will look for such effects with regard to 
the QPCs. 



C. Simulation of the QPC 

In this section we will first obtain the bare confinement 
potential created by the QPC for the geometry given in 
figure^ and then go on to discuss the effects of screening. 
The potential generated by such gates can be calculated 
by the scheme proposed by Davies et alm^. 

The model parameters are taken from the relevant ex- 
perimental samples2il4, where the applied gate voltage is 
—0.3V, the width at the tip is about 200nm, and the tip 
separation Ay ^ 300nm. The 2DES is taken to be 85nm 
below the surface. 

We define the QPC using rectangles and polygons 



which are shown in figure ^ as red (dark) and white ar- 
eas. In figure 13 we show the bare confining potential for 
the parameters given above. The electrons are filled up 
to the Fermi energy [Ep ~ 7meV, corresponding to a 
typical electron surface density Uei ~ 1.7 x 10^^ cm^^). 
Using such parameters, the full screening calculation to 
be discussed below will reveal the electrons to be depleted 
beneath the QPC, say at all the dark (blue) regions in 
FigH 

In our numerical simulations, we have mapped the 
unit cell containing the QPC of physical dimensions 
3.3^mxl.8/im to a matrix of 200 by 200 mesh points 
in the absence of a magnetic field and l.l/xmxl.8^m to 
a matrix of 48 by 96 mesh points in the presence, which 
allows us to perform numerical simulations within a rea- 
sonable computation time. With regard to numerical ac- 
curacy, we estimate that, for typical electron densities, 
the mean electron distance, i.e. the Fermi wavelength, 
is larger than 40nm. Hence, the number of mesh points 
considered here allows us to calculate the electron den- 
sity with a good numerical accuracy. We also performed 
calculations for finer meshes and the results do not dif- 
fer quantitatively (at the accuracy of line thicknesses), 
whereas the computational time grows like the square 
of the number of the mesh points. We should also note 
that due to computation time concerns we had to use a 
smaller unit cell in the presence of the magnetic field, 
which yields finite size effects close to the boundaries of 
the sample (e.g. see Fig[7l3). The features observed are, 
in principle, negligible and they tend to disappear when 
the unit cell is taken to be larger and mapped on a larger 
matrix. 

We now discuss the resulting bare and screened poten- 
tial for a realistic QPC defined by surface gates, with a 
tip opening Ay = 300nm. Figure |21 represents the exter- 
nal potential created by the QPC gate structure at the 
surface, calculated in the plane of the 2DES located at 
z = 85nm below the surface, with an applied potential 
—0.3V. In the upper panel we show a 3D plot and a pla- 
nar projection, together with four guide lines, which in- 
dicate the location of the cross-sections through the that 
are displayed in the lower panel. The level of the Fermi 
energy of the system (to be assumed below) is indicated 
in the 3D plot as well. These results have been obtained 
numerically from Eq. The barrier is formed by the 
regions of elevated potential. 

At the first glance one observes that the potential land- 
scape is smoothly varying. This is purely an effect of the 
relatively large distance to the gate, as screening effects 
have not yet been included. For the given Fermi energy 
(obtained from the electron density in the bulk) and the 
tip separation. Ay > 100 nm, the number density of elec- 
trons inside the QPC opening fulfills the validity relation 
of the TFA, i.e. ne;(center)aB ^ 1. At the positions 
where the height of the barrier becomes larger than the 
Fermi energy (light line in the 3D plot and horizontal 
dashed line in the lower panel), the probability to find 
an electron is zero within the TFA. 
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FIG. 3: The bare confinement potential generated by the 
QPC, defined by the polygons shown in figure The color 
scale indicates the strength of the confinement. 



We proceed in our discussion with a comparison of the 
screened potential shown in Fig. 21 to the bare confine- 
ment potential discussed up to now (Fig. 13). The self- 
consistent potential is obtained from the formalism de- 
scribed above for periodic boundary conditions at zero 
temperature and zero magnetic field. The electrons arc 
filled up to the Fermi energy (shown by the gray thick 



line on the surface of the color plot and dashed line in 
the lower panel), such that no electrons can penetrate 
classically into the barrier above those lines. The first 
observation is that the potential profile becomes sharper 
for the screened case and strong variations are observed 
in the vicinity of the QPC. These shouldcr-likc local max- 
ima near the QPC represent the same feature seen in the 
example of the square barrier discussed previously, and 
we have pointed out that they stem from g-dependent, 
non-linear screening. This will become more important 
when we consider a magnetic field, since the local "pin- 
ning" of the Landau levels to the Fermi energy in these 
regions will produce compressible regions surrounded by 
incompressible regions. 

An interesting feature occurs near the opening of the 
QPC, namely a local minimum which is a result of the 
non-linear screening. We point out that somewhat simi- 
lar physics has been found (using spin-dcnsity-functional 
theory^) to lead to the formation of a local bound state 
inside a QPC, which has been related to the "0.7" 
anomaly, linking it with Kondo physics. We believe this 
feature to be a very important result of the self-consistent 
screening calculation, and we note that it may affect 
strongly the transport properties of the QPC both in the 
presence or absence of a magnetic field. We will discuss 
the influence of this local minimum on the formation of 
the incompressible strips in section llV Dl where wc calcu- 
late the density and potential profiles including a strong 
perpendicular magnetic field. 

It is known from the experiments that the interference 
pattern and the transmission properties strongly depend 
on the structure of the QPCs, such as the distance of 
the 2DES from the surface, the applied gate voltage, the 
sharpness and the geometry of the edges, as well as the 
width of the opening of the QPC. The effect of the first 
two parameters can be understood by following the sim- 
ple arguments of linear screening as shown for the square 
gate model: if the distance from the QPC to the 2DES 
increases, the potential profile becomes more and more 
smooth. The screened potential changes linearly with 
the applied gate potential (see Ea. lfTT|) '). The geometric 
parameters have to be adapted to the experiment in ques- 
tion. Note that the shape of the QPCs has already been 
discussed in the literature (see Ref. |2 and references con- 
tained therein). The effect of the size of the QPCs, how- 
ever, has not been considered for large Ay (> lOOnm), 
and we believe this to be an important parameter for the 
interferometer experiments. 

We start our investigation by looking at the opening 
of the QPC with increasing tip separation of the metal 
gates used to define the QPC. In this section, we work 
at zero temperature and magnetic field, with a constant 
bulk electron density. 

In figureElwe depict the self-consistent potential at the 
center of the QPC (j/ = 550nm), while changing the tip 
separation (Ay) between 100 and 500nm. We see that 
for the narrowest separation the potential profile looks 
rather smooth and a minimum is observed at the center. 
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FIG. 5: The screened potential at x = 550nm for five different 
tip distances. Note that the x interval used for the calculation 
smaller than in the previous figure, since we concentrate on 
the bulk structures rather than the edge ones. 
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FIG. 4: The screened potential (upper panel) seen by a 2DES 
at 85nm below the surface and some characteristic cuts along 
the X- axis, together with an indication of the Fermi level Ef 
(lower panel). The color scale represents the strength of the 
potential, and the cross-sections are indicated by the same 
line code as in Fig. |3 



If we increase Ay(^ 300nm) we see that the screening be- 
comes stronger, leading to more pronounced shoulders on 
the sides and a deeper minimum at the center. For even 
larger separations (Ay > 300nm) a local maximum starts 
to develop at the center, since the electrostatic potential 
energy is no longer strong enough to repel the electrons 
from this region. Basically all the non-linear features 
observed result from the competition between the gate 
potential, which simply repels the electrons, and the mu- 
tual Coulomb interaction, i.e. the Hartree potential. It 
is obvious that for narrower tip separations only a few 
electrons will remain inside the QPC opening and there- 
fore TFA type approximations will not be justified any 
longer. 

Summarizing this section, we have determined the 
screened potential profile in a realistic QPC geometry, 
pointing out features resulting from non-linear screen- 
ing. We have observed that a local extremum occurs 
at the center of the QPC, and have traced the depen- 
dence on the width Ay between the QPC tips. These 
features, as mentioned before, become more interesting 
if a magnetic field is also taken into account, where they 
lead to stronger spatial inhomogeneities in the electron 
distribution. Our next step is thus to include a strong 
quantizing perpendicular magnetic field and examine the 
distribution of the incompressible strips where the im- 
posed external current is confined^SiS^. 



D. Finite temperature and Magnetic field 



Once the initial values of the screened potential and the 
electron distribution have been obtained for T = 0, B = 
0, using the scheme described above, one can calculate 
these quantities for finite field and temperature as fol- 
lows: replace the zero temperature Fermi function with 
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FIG. 6: [a]-[d] Color-coded plot of the local filling factor ver- 
sus position {x, y) for a square sample of width = ay = 
0.8 /xm; white indicates ^{x, y) — 2. The average density is 
taken to be 3.0 ■ lO" cm"^; kBT/Ep = 0.02. [e] A sketch of 
the Hall resistance as a function of magnetic field 



the finite temperature one and insert the bare Landau 
DOS 



Db{E) = —Y^5{E~ En), En = hLo,{n+l/2) (12) 

n=0 

into equation Q instead of Dq. In our numerical scheme 
we first start with relatively high temperatures (i.e. a 
smooth Fermi function) and then decrease the tempera- 
ture slowly until the desired temperature is reached. A 
Newton- Raphson method is used for the iteration process 
and at every iteration step the electro-chemical potential 
is checked to be constant. 

Before proceeding with the investigation of the QPC 
geometry at B > 0, we would like to make clear the re- 
lation between the quantum Hall plateaus and the exis- 
tence of the incompressible strips following the arguments 
of Siddiki and GerhardtsS^. Fig. presents the local fill- 
ing factors of a relatively small Hall bar. together with an 
illustrative Hall resistance curve. At the high magnetic 
field side ((b), iy{0,0) < 2) there are no incompressible 



strips, thus the system is out of the Hall plateau. When 
approaching from the high B side to the plateau a single 
incompressible strip at the center develops. When the 
width of this strip becomes larger than the Fermi wave- 
length, the system is in the quantum Hall state ((d), 
1/(0, 0) = 2). If we decrease the field strength further the 
center incompressible strip splits into two and moves to- 
ward the edges ((c) !^(0,0) > 2). As long as the widths 
of these strips are larger or comparable with the Fermi 
wavelength the system remains in the plateau. This is 
the regime in which an interferometer may be realized. 
Further decreasing the magnetic field leads to narrower 
incompressible strips which finally disappear if the widths 
of them become smaller than the average electron dis- 
tance. Then the system leaves the quantized plateau. 
The distribution of the incompressible strips and the on- 
set of the plateaus, of course, depends on the disorder 
potentials^ and the physical size of the sample. However, 
the experiments considered here are done using narrow 
and high mobility structures, thus the above scheme will 
cover the experimental parameters. 

In this subsection we present some of our results ob- 
tained within the TFA using periodic boundary condi- 
tions, considering two different tip separations, while 
sweeping the magnetic field. First we will fix the gate 
potential to Vqpc = —0.3V and sweep the magnetic field 
for Ay = lOOnm, while keeping the electron number den- 
sity, i.e. the Fermi energy, constant. Second we examine 
the potential profile for Ay = 300nm and comment on 
the possible effects on the coherent transport properties. 

In figure[71wc plot the local filling factor (i.e. the nor- 
malized density) distribution of the 2DES projected on 
the ccy-plane, together with the same quantity for some 
selected values of y, at average filling factor {D) one. 
From the y = Onm curve (solid lines) in FigEIb], one 
can see that the electrons beneath the QPC arc depleted 
(shaded, dark (blue) regions) (300 < x < 800nm), while 
the electron density reaches finite values while approach- 
ing the opening of the QPC {y ~ 850nm). At D ~ 1 
one does not observe any incompressible regions, since 
the Fermi energy is pinned to the lowest Landau level. 
Hence the electron distribution is rather smooth and the 
current distribution will just be proportional to the num- 
ber of electrons, similar to the Drude approach. For this 
case the external potential is screened almost perfectly 
and the self-consistent potential is almost flat, thus one 
can assume that the corresponding local wave functions 
are very similar to the ground state Landau wave func- 
tions. 

The first incompressible region occurs when the Fermi 
energy falls in the gap between two low-lying Landau 
levels. Then the electrons exhibit a constant density 
and thus cannot screen the external potential. In fig- 
ure ©a, we show the electron distribution for D = 1.1. 
The black regions denote a local density corresponding 
to filling factor i/ = 2, which does not percolate from 
the left side of the sample (which we might identify with 
the source) to the right side (drain). Here one can see 
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FIG. 7: [a] The top view of the local filling factor, v{x,y), 
distribution of the 2DES, for average filling factor one, in the 
plane located at 2 = 85nm below the surface, at the "de- 
fault" temperature, kT/huJc = 1/50, which will be used in 
all subsequent plot. The color scale depicts the density of 
electrons, whereas the dark shaded areas indicate the elec- 
tron depleted regions, [b] Side view of the local filling factor 
for y = Onm(solid line), 468nm (dashed line), 660nm (dot- 
ted line), 750nm (dash-dotted line) and 900nm (dash-dotted- 
dashed line) . The horizontal lines in [a] shows the positions 
of the cuts in [b], with the same line code. Note that the den- 
sity has local minima at large and small x, which are finite 
size efi'ects mentioned in the text. 



well developed incompressible puddles, at the regions 
150nm< X < 250nm, Onm< y < 450nm (and four other 
symmetric ones) and two smaller puddles at the entrance 
of the QPC. These structures will remain unchanged even 
if one considers a larger unit cell, since they manifest the 
q- dependency, i.e. the rapid oscillations of the Fourier 
transform of the confining potential of the QPC. 

In these regions the self-consistent potential exhibits 



a finite slope. Accordingly the wave functions will be 
shifted and squeezed, i.e. they are now superpositions 
of a few high order Landau wave functions with re- 
normalized center coordinates. This behavior has been 
showni2Si2& for the translationally invariant model. Here 
we did not include the finite extent of the wave functions, 
here to avoid lengthy numerical calculations. 

The incompressible regions shift their positions on the 
xy plane depending on the strength and the profile of 
the confining potential. In figure JHJb we show the filling 
factor distribution where the bulk filling factor is almost 
two. We see that four incompressible strips arc formed 
near the QPC. However the QPC opening remains in 
a compressible state, with local filling factor less than 
two, where we expect that the self-consistent potential 
is essentially flat. Further increasing the average filling 
factor, we observe that the bulk becomes completely com- 
pressible and two incompressible strips are formed near 
the QPC which percolate from bottom to top, creating a 
potential barrier with a height of hujc, see Fig. (jS);). For 
even higher filling factors, they merge at the center of 
the QPC (Fig.|Hli). In that case, the potential within the 
QPC will then no longer be flat, due to poor screening. 
We should also note that for a small width Ay of the 
QPC opening, merging of the incompressible strips will 
happen only in a very narrow B interval, and a quan- 
titative evaluation within our TFA can not be always 
satisfactory, as the number of electrons inside the QPC 
becomes too low. Further decreasing the field strength 
(increasing the average filling factor) results in two sep- 
arate incompressible strips winding around the opposite 
gates making up the QPC, as shown in Fig. Thus, 
dissipationless transport through the QPC, with a quan- 
tized conductance, becomes possible. At the lowest field 
values considered in this figure, we see that the innermost 
incompressible strips (with 1/ — 2) become smaller than 
the Fermi wavelength and thus they essentially disappear 
and no longer affect the transport properties. This point 
has been discussed in detail in a recent work by Siddiki 
et alm^. The scheme described above now starts to re- 
peat, but with incompressible strips having a local filling 
factor of 4. 

We now discuss the effects of increasing the separation 
parameter, which we choose to be Ay = 300nm in figure 
© . At the strongest magnetic field ^j^) , only very small 
regions are incompressible and the electron distribution is 
similar to Fig|SJi, where the incompressible regions result 
from local unpinning of the Fermi energy from the low- 
est Landau level due to dependent screening, i.c the 
shoulder-like variation of the potential near the QPC dis- 
cussed earlier. By decreasing B, an interesting structure 
is observed at the center of the QPC: an incompress- 
ible island. In figure ^jp, we have tuned the magnetic 
field such that the bulk of the 2DES is incompressible, 
meanwhile the entrance to the QPC remains compress- 
ible. The strong variation of the self-consistent potential 
at the center of the QPC can generate a pronounced effect 
on the current passing through the QPC (see figure [TUI 
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FIG. 8: The local filling factor distribution for different av- 
erage filling factors {v), which is defined by the number of 
the electrons in the unit cell, [a] i> = 1.1 [b] 1.2 [c] 1.4 [d] 
1.6 [e] 1.8 [f] 2.2. The color scale depicts the local electron 
concentration, whereas the abrupt colors indicate the even- 
integer filling factors, i.e. incompressible strips, (black for 
v{x,y) = 2, magenta for i/{x,y) = 4. 



and the related text). For a lower magnetic field strength 
the distribution of the incompressible region is just the 
opposite (c). Now we see a large compressible puddle at 
the center, surrounded by incompressible regions, which 
can percolate from source to drain. Coherent, dissipa- 
tionless transport can be expected in this case. Further 
decreasing the magnetic field we observe that the struc- 
ture is smeared out and the tip region becomes com- 
pressible, nevertheless there are two large incompressible 
regions close to the entrance of the opening. The two 
incompressible strips wind around the gates, as shown 
in Fig. Ol. Finally, a scheme similar to that observed 
earlier in figure |Hli-e is also seen now, for Ay = 300nm. 

Another remark which we would like to make concerns 
the edge profile of the sample itself and of the QPC. It has 
been shown both cxperimentallji^ and thcorcticalljiSSi^ 
that for an (almost) infinite potential barrier at the edges 
of the sample, the Chklovskii edge state picture 
breaks down, i.e. no incompressible strips near the edge 
can be observed. Meanwhile for smoothly varying edge 
potential profiles many incompressible strips are present, 
if the bulk filling factor is larger than two (for spinless 
electrons: four). We believe that, within the MZI setup 
both of these edge potential profiles might co-exist. At 
the edge regions of the sample, where lateral confine- 
ment is defined by physical etching, the potential profile 
differs from of the one generated by the top gates, due 
to different separation thicknesses and also lateral sur- 
face charges generated by etching. In principle gate and 
etching defined edges impose different boundary condi- 
tions, and the effects on screening at a 2DES have been 
discussed before24. These two profiles will certainly af- 



fect the group velocity, since the slope of the potential 
depends on the (lateral) boundary conditions. Follow- 
ing the arguments of Ref . [2^l25| . which essentially pre- 
dict that the dissipative current is confined to the in- 
compressible strips, the widths of these strips will also 
define the slope, hence the velocity of the electrons will 
be determined by the edge profile. The velocity of the 
edge electrons were investigated experimentallji^ and the 
magnetic field dependency was reported as B^^^. There 
it was noted that a self-consistent treatment is necessary 
to understand their findings, which we would like to dis- 
cuss in a future publication. 

The important features to note in these results are (i) 
in general, electron-electron interactions have a remark- 
able effect, leading to the formation of a local extremum 
in the potential at the center of the QPC, which even at 
low electron densities seems to be well described by the 
TFA; (ii) the narrow compressible/incompressible strips 
formed near the QPC are a direct consequence of the 

dependent screening. 



E. Comments on coherent transport 

A complete calculation of coherent transport requires 
a deeper analysis of the wave functions and is beyond 
the scope of this work, which has been devoted to self- 
consistent realistic calculations of the potential and den- 
sity profiles. In principle, one can follow the arguments of 
the well developed recursive Green's function technique^ 
in the absence of magnetic field and the method devel- 
oped recently even in the presence of a strong field^. 

Instead we would like to examine the potential dis- 
tribution across the QPC and comment on the possible 
effects of interaction on the wave functions, and thereby 
(indirectly) on transport. In figure [TUI we depict the po- 
tential profile across the QPC for the parameters used to 
obtain figure |S1 As expected for D = 1.0 (dashed (red) 
line) the 2DES is "quasi" metallic, hence the external 
potential is perfectly screened, and the wave functions 
are left almost unchanged. The two incompressible is- 
lands seen at the entrance of the QPC in Fig|Hli lead to a 
minor variation of the screened potential at a; = 300nm 
and X = 800nm, depicted by the solid (black) line for 
D ~ l.l. A drastic change is observed when the bulk be- 
comes incompressible {D = 1.2) and the opening remains 
compressible: Now the 2DES cannot screen the exter- 
nal potential near the openings of the QPC, where we 
see a strong variation. The strong perpendicular mag- 
netic field changes the potential profile near the QPC 
via forming incompressible strips, and local minima are 
observed at the entrance and the exit. In these regions 
the electrons are strongly localized and the wave func- 
tions are squeezed. The situation is rather the opposite 
for D ~ 1.4, where two incompressible strips located near 
the QPC, formed due to g-depcndent screening, merge 
at the opening. One observes a barrier with the height 
of fkUc, which essentially is a direct consequence of the 
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FIG. 9: The local filling factor distribution for different av- 
erage filling factors (i?), for a tip separation 6y = 300 nm. 
Note that, the number of electrons in the unit cell is changed, 
since the depleted areas are larger than of the previous fig- 
ure, [a] v = 1.14 [b] 1.2 [c] 1.34 [d] 1.4, [e] 1.6 and [g] 
3.1. The color scale depicts the local electron concentration, 
whereas the high-contrast color regions indicate the even- 
integer filling factors, i.e. incompressible strips, (black for 
v{x,y) — 2, magenta for v{x,y) = 4). The calculations 
are done at ksT/huJc = 1/50 for an average electron density 
1.7 X 10""cm-^ 



incompressible strip at the center and electrons have to 
overcome this barrier. Further decreasing the magnetic 
field smears out the barrier gradually, until the system 
becomes completely compressible and wc arc back in the 
case of figurelTUh (also with regard to the transport prop- 
erties). 



V. SUMMARY 

The study that motivated the present authors 
was the Quantum Hall efi^ect based Mach-Zender 
interferometesSii. The results of these experiments de- 
viate markedly from the naive single-particle Landauer- 
Biittiker picture, and interaction effects have been sug- 
gested as a possible explanation. These effects may in- 
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FIG. 10: The self-consistent potential profile across the QPC, 
plotted for characteristic values of the average filling factor. 
Calculations are done at the default temperature and electron 
density. 



elude decoherence due to potential fluctuations brought 
about by electron-electron or electron-phonon interac- 
tions (together with other noise sources). A more de- 
tailed understanding of electron-electron interactions in 
this setup, as well as of those features of the interferom- 
eter that are specific to the physics of the Quantum Hall 
effect, hinges on an analysis of the self-consistent static 
potential landscape near the QPCs, which represent the 
most crucial components of the setup. 

Therefore, in this work, we have taken into account 
the electron-electron interaction within the TFA, consid- 
ering realistic geometries of QPCs, calculating the self- 
consistent potential and electronic density profiles, and 
commenting on possible effects on transport. 

The outcome of our model calculations can be sum- 
marized as follows: (i) We have obtained the electro- 
static potential profile for the QPC geometries used in 
the experiments by solving the Laplace equation semi- 
analytically. (ii) We have demonstrated for a simple 
square well barrier that the screened potential in a 2DES, 
even in the absence of a magnetic field, is strongly depen- 
dent on the initial potential profile and on the distance 
between gates and 2DES. (iii) The screened potential has 
been calculated within the TFA for a QPC at vanishing 
field, where we have observed two interesting features: 
a local extremum at the center of the QPC and strong 
shoulder-like variations near the QPC. (iv) In the pres- 
ence of a magnetic field, the formation and the evolu- 
tion of the incompressible regions has been examined and 
three different cases have been observed: (a) the system 
is completely compressible, (b) An incompressible re- 
gion and/or strip, which does not percolate from source 
to drain, generates a local extremum at the entrance/exit 
of the QPC (c) The center of the QPC becomes incom- 
pressible, with or without a compressible island, hence 
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the incompressible strip percolates from source to drain. 
We note that the local minimum found at the center of 
the QPC for certain tip separations, being a clear inter- 
action effect, coincides with the findings of Hirose et al 
0. 
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